Deeply Virtual Compton Scattering Beam-Spin Asymmetries 
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The beam spin asymmetries in the hard exclusive electroproduction of photons on the proton 
(ep —> ep7) were measured over a wide kinematic range and with high statistical accuracy. These 
asymmetries result from the interference of the Bethe-Heitler process and of deeply virtual Compton 
scattering. Over the whole kinematic range {XB from 0.11 to 0.58, from 1 to 4.8 GeV^ -t from 
0.09 to 1.8 GeV^), the azimuthal dependence of the asymmetries is compatible with expectations 
from leading-twist dominance, A ~ a sin 0/(1 + c cos 0). This extensive set of data can thus be used 
to constrain significantly the generalized parton distributions of the nucleon in the valence quark 
sector. 

PACS numbers: 12.40.Vv,13.40.Gp,13.60.Fz,13.60.Hb,13.60.-r,14.20.Dh,24.85.+p 



The structure of the nucleon, the lightest of all bary- 
onic states, has been studied in the past using two com- 
plementary approaches. Elastic electron scattering mea- 
sures form factors which reflect the spatial shape of 
charge distributions while deep inelastic scattering 
provides access to parton distribution functions that en- 
code, in a fast moving nucleon, the momentum fraction 
carried by the constituents Q]- The formalism of Gener- 
alized Parton Distributions (GPDs) 0, 0, 0] unifies these 
approaches and provides much greater insight into nu- 
cleon structure @, Q, through the coherence between 
states of different longitudinal momentum fractions, the 
correlation between transverse coordinates and longitu- 
dinal momentum of the partons Q, the distribution of 
forces exerted upon partons P| (information inconceiv- 
able to obtain just a few years ago) and the angular mo- 
mentum carried by each type of parton 

Deeply virtual Compton scattering (DVCS) on the pro- 
ton (7*p — > 7p), in the Bjorken regime where the pho- 
ton scattering occurs at the quark level, is the process of 
choice to attain an experimental determination of GPDs. 
Pioneering observations of DVCS [13, [ll|, [H, [H, Ei, EE , 
16[ , though of limited experimental accuracy, are all com- 
patible with a description of the observables in terms of 
GPDs, both in the gluon and in t he q uark sector. More- 
over, a recent precise experiment [17| gave good indica- 
tions of the onset of scaling in this process at relatively 
modest values of the 7* virtuality. 

In this context, this work presents the first system- 
atic and precise exploration of a sensitive observable, the 
beam-spin asymmetry of the reaction ep — > ep"f. Neglect- 



ing a twist-3 DVCS term, this asymmetry arises from the 
interference between the Bethe-Heitler (BH) and DVCS 
processes (that is, where the photon is emitted by the 
electron or by the target nucleon, respectively). At lead- 
ing twist, it is primarily sensitive to the imaginary part of 
the DVCS amplitude and thus to a specific linear combi- 
nation of the proton GPDs H, H and E, with arguments 
x{= ±f), ^ and t. Each proton GPD involves a weighted 
sum over the quark flavors. The beam-spin asymmetry 
is defined as 



(1) 



where the arrows correspond to beam helicity -1-1 and —1. 
It depends on Q^, xb, t, defined in Fig. [1] and on the 
angle 4> between the leptonic and hadronic planes. Har- 
monic decompositions of the cross sections d^cr, divided 
among contributions from BH, DVCS and interference 
(INT) terms, have been proposed 18, 19]. In the nota- 
tion of Ref. Ilii, the cross sections, up to some kinematic 
factors, can be expressed in terms of the (/)-harmonics 
cos n(f> and sinn0, with n from to 3 and S — BH, 
INT or DVCS. At the twist-2 level, which according to 
Ref. |l7[ is largely dominant at least up to |t| = 0.35 



GeV^, the numerator of Eq. ([T]) gets a contribution from 
only, while the denominator contains the coefficients 
and c^^°^, in addition to c™(n = 0,1,2) 
calculable in QED in terms of the proton elastic form 
factors. At leading twist, one obtains 
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FIG. 1: (Color online) Schematic representation of the 
leading-order handbag diagram contribution to DVCS, where 
X is the average longitudinal momentum fraction of the ac- 
tive quark in the initial and final states (measured in terms 
of the average hadron momentum (p -\-p )/2), while 2^ is 
their difference; it is related to the Bjorken scaling variable 
by^ ~ xb/{2 — xb)- The squared four-momentum transfer to 
the target is t = {p' — p)^ , and the squared four-momentum 
of the virtual photon is —Q^. 



where the parameters a, c and d may be expressed in 
terms of the above mentioned harmonic coefficients. The 
DVCS and INT harmonic coefficients may in turn be 
written in terms of Compton form factors related to the 
corresponding GPD by 



2x 



^2 -a;2 



ff(x,e,i) , (3) 



-SmW = Hi^,tt)-H{-C,^,t) , (4) 

TT 

up to corrections of order of the strong coupUng constant, 
with similar expressions for H, £ and £. The GPD H 
yields the dominant contribution to the harmonic coef- 
ficients considered above. Neglecting the small contri- 
butions from the three other GPDs, one can express the 
beam-spin asymmetry A in terms of only 3?e7i and QuiH. 
Thus in this approximation, which is expected to hold for 
small values of |t | , the parameters a, c and d of Eq. ([2]) are 
uniquely related to the imaginary and real parts of the 
Compton form factor H, yielding respectively the GPD 
H at points x = ±^ and the principal value integral of 
Eq. ([3]). Going beyond this approximation requires ad- 
ditional theoretical or experimental constraints on the 
other GPDs. 

The experiment took place in Hall B of Jefferson Lab, 
using the CEBAF 5.77 GeV electron beam (with average 
polarization P = 0.794), a 2.5 cm-long liquid-hydrogen 
target and the CLAS spectrometer The three final- 
state particles from the reaction ep ejry were detected. 
For this purpose, a new inner calorimeter (IC) was added 
to the standard CLAS configuration, 55 cm downstream 
from the target, in order to detect 1 to 5 GeV photons 
emitted between 4.5° and 15° with respect to the beam 
direction. This calorimeter was built of 424 tapered lead- 
tungstate crystals, 16 cm long and with an average cross- 
sectional area of 2.1 cm^, read out with avalanche photo- 
diodes and associated low- noise preamplifiers. The whole 



IC was operating at a stabilized temperature of 17°C, and 
monitored with laser light homogeneously distributed on 
all crystals. The calorimeter was calibrated several times 
during the run using the two-photon decay of neutral pi- 
ons. Energy and angle resolutions of 4.5% and 4 mrad 
(for 1 GeV photons) were achieved. In conjunction, a 
specifically-designed superconducting solenoid was used 
to trap around the beam axis the background originating 
from M0ller electrons, while permitting detection of the 
recoil protons up to 60° . 

Events were selected if an electron had generated the 
trigger, one and only one proton was identified and only 
one photon (above an energy threshold of 150 MeV) was 
detected in either the IC or the standard CLAS calorime- 
ter EC. Electrons were identified through signals in the 
EC and in the Cerenkov counters. From time-of-flight 
information, track length and momentum, protons were 
unambiguously distinguished from positive pions over the 
whole momentum range of interest. All clusters detected 
in the IC were assumed to originate from photons, while 
additional time-of-flight information was used in the EC 
to separate photons from neutrons. For all three final- 
state particles, fiducial cuts were applied to exclude de- 
tector edges. 

Operating at a luminosity of 2 x lO"^** cm^^g"! (a record 
for CLAS), the accidental coincidences were negligibly 
small, as well as the pile-up probability in the IC, ex- 
cept for the most forward photons below 6°. Events 
considered here include the kinematic requirements : 
> 1 GeV^, 7*p invariant mass W > 2 GeV and 
scattered electron energy E' > 0.8 GeV. The mere se- 
lection of the three final-state particles results in the ob- 
servation of characteristic peaks in distributions of all 
kinematic variables expressing the conservation of total 
four-momentum in the reaction ep — *■ epj, as exemplified 
by the dotted curves in Fig. [2j Requiring in addition a 
missing transverse momentum smaller than 0.09 GeV, an 
angle between the 7*p' and jp' planes smaller than 1.5°, 
a photon detected within 1.2° of the direction inferred 
from the detected electron and proton, and a maximal 
missing energy Ex of 0.3 GeV, results in clean peaks for 
the events of interest. These kinematic cuts are to some 
extent redundant (except for the background to be dis- 
cussed below) and are quoted here for the case where the 
emitted photon is detected in the IC, that is for 92% of 
the events. In the case of photons detected in the EC, 
these cuts are about twice as large because of the poorer 
resolution. 

In spite of this selection, a contamination of events 
originating from the ep epw^ reaction, followed by the 
subsequent asymmetric decay of the neutral pion, is al- 
ways possible. For these events, one of the photons is 
not detected, because it is either below threshold or out- 
side the calorimeters' acceptance. This physical back- 
ground is estimated using the number iV^o of measured 
ep — > epir'^ events, identified unambiguously when the 
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FIG. 2: (Color online) Distributions in cone angle O-^y for the 
ep — > epY reaction (left) and in missing energy Ex for the 
ep —> ep'yX reaction (right), before (black dotted curve) and 
after (red solid) all kinematic cuts discussed in the text but 
the one on the histogrammed variable, given by the location 
of the arrow. The thin solid black line represents the physi- 
cal background, calculated from measured ep epn^ events. 
The distributions are integrated over all kinematic variables 
and apply to the case where the photon is detected in the IC. 



two photons are detected [2l| , and multiplying it by the 
ratio of acceptances Acc^o/Acc^o, where the "I7" ac- 
ceptance is to be understood with the photon satisfy- 
ing all the ep — > epj event selection cuts. This ratio, 
which depends mostly on the photon geometrical cuts 
and on the relevant resolutions, has been calculated with 
the standard CLAS simulation package and a simplified 
fast Monte-Carlo, the two results being used to evaluate 
the corresponding systematic uncertainties. The back- 
ground proportion / varies between 1 and 25% depend- 
ing on the kinematic bin, 5% in average. The number of 
ep — > epj events is then, for each beam-helicity state and 
for each elementary bin in the four kinematic variables 
(see below), N = Nep^epiX - (Acc^'2/Acc^^)7V^?, and 
the asymmetry A ^ {N - 'N)/P{N + "N). Finally, radia- 
tive corrections were applied [22] . These tend to increase 
the asymmetries very slightly. 

The data were divided into thirteen bins in the {xb, 
Q^) space as per Fig. [31 five bins in — t (defined by 
the bin limits 0.09, 0.2, 0.4, 0.6, 1 and 1.8 GeV^) and 
twelve 30° bins in 0. Bin-size corrections were applied. 
Whether integrated in t or in each i-bin (Fig. [3]), the 0- 
distributions were always found to be compatible with 
Eq. ^ with d = 0. The parameter d is expected to 
be smaller than 0.05 over our kinematic range, and in- 
deed was found compatible with zero, within statistical 
accuracy, when including it in the fit. The deviation 
from a pure sine function as \t\ increases is seen in all 
{xb,Q'^) bins and results in the parameter c becoming 
negative [2^. The parameter a is the best estimate of 
A{90°) and is represented in Fig. [H Point-to-point sys- 
tematic uncertainties arise mostly from the background 
subtraction: AAt, — {A — A.^o)Af /{I — /), where the 
relative error on / is conservatively estimated to be 30% 
and A^o is the asymmetry for the reaction ep —> epir'^, 
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FIG. 3: (Color online) Left: kinematic coverage and binning 
in the (x_b, Q^) space. Right: ^(0) for 2 of the 62 (a;_B, Q^, 
t) bins, corresponding to {xb) = 0.249, (Q^) = 1.95 GeV^ 
and two values of (t). The red long-dashed curves correspond 
to fits with Eq. ([2]) (with d = 0). The black dashed curves 
correspond to a Regge calculation The blue curves corre- 
spond to the GPD calculation described in the text, at twist-2 
(solid) and twist-3 (dot-dashed) levels, with H contribution 
only. 



ranging between 0.04 and 0.11 at 90° [21|. The sensi- 
tivity of the results to the event selection cuts was stud- 
ied as well. From these two sources of information, the 
systematic uncertainty on a was inferred to be 0.010, in- 
dependent oi Xb, and t. An overall normalization 
uncertainty arises from the uncertainty in the beam po- 
larization (3.5%). Additional details on the experiment 
and on the data analysis may be found in Ref. 2^. 



The wide kinematic coverage of the present data is 
important for global analyses of ep — > ep^ observables 
and for a model-independent extraction of DVCS ampli- 
tudes. The beam-spin asymmetries are especially, but 
not uniquely, sensitive to the GPD H . When combined 
with other observables more sensitive to H and E, as 
well as with unpolarized cross sections, it will be possi- 
ble to obtain the real and imaginary parts of the Comp- 
ton form factors of all GPDs, as defined in Eqs. ([3]) 
and Additional theoretical work is also required, to 
clarify how power-suppressed contributions not included 
in Ref. ^] would affect the relations between observables 
and GPDs [25[ . Presently, GPDs may be calculated using 
theoretical models based on constituent quarks, on a chi- 
ral quark-soliton description of the nucleon, on light-cone 
or other frameworks. The first moments of GPDs are be- 
ing calculated using lattice QCD techniques. But none of 
these calculations are developed to the point of making 
the link to DVCS observables. Alternatively, constrained 
parameterizations have been used to make predictions of 
DVCS beam-spin asymmetries. Following Refs. 2^ 27 1, 
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FIG. 4: (Color online) a = y4(90°) as a function of -t. Each 
individual plot corresponds to a bin in {xb,Q^)- Systematic 
uncertainties and bin limits are illustrated by the grey band in 
the lower left plot. Black circles are from this work. Previous 
results are from Ref. [121 (red square) or extracted from cross 
section measurements [17| (green triangles), at similar - but 
not equal - values of {xb) and (Q^). See Fig. [3] caption for 
curve legend. 



such a parameterization of the GPD H may be 



1 .1-1/31 

d/3 / daS{x-(3- ^a)h'' (/3, a, t) 

2 ^ 

(5) 



with a, t) = qiP)nbif3, a)e~°''^^^-'^^' 



(6) 



where and q{(3) are the electric charge and unpolar- 
ized parton distribution for quark flavor g, tt;, a profile 
function and a[ is a Regge slope adjusted to recover 
the proton form factor Fi from the first moment of the 
GPD. Eq. (HI) extends the ansatz of Ref. |23| for the t de- 
pendance to non-zero values of The D term in Eq. ([5]) 
is calculated within a quark-soliton chiral model Q • Us- 
ing predetermined parameters, the calculations of beam- 
spin asymmetries yield the solid and dot-dashed curves 
in Figs. [3]and|4l without and with a twist-3 term calcu- 
lated in the Wandzura-Wilczek approximation Q. The 
predictions overestimate the asymmetries at low \t\, es- 
pecially for small values of xb and/or Q^. Variations 
of the parameter b entering the profile function TTf, do 
not resolve this problem, which may indicate that dou- 
ble distributions are not flexible enough to reproduce this 
behaviour. 



Alternatively, description of the process in terms of 
meson (or more generally Regge trajectory) exchanges 
has been attempted (j^. [28|. DVCS may be viewed as p 
production followed by p-j coupling in vacuum or in the 
nucleon field. In addition to pole contributions in the t 
channel [2^ , the box diagram that takes into account p- 
nucleon intermediate states has been evaluated [l^l . This 
calculation, represented by the dashed curves in Figs. [3] 
andlH is in fair agreement with our results up to = 2.3 
GeV^. The significance of this dual description (Regge 
vs. handbag) remains to be fully investigated. 

In summary, the most extensive set of DVCS data to 
date has been obtained with the CLAS spectrometer, 
augmented with specially designed small-angle photon 
calorimeter and solenoid. Beam-spin asymmetries were 
extracted in the valence quark region, as a function of all 
variables describing the reaction. Present parameteriza- 
tions of GPDs describe reasonably well, but not perfectly, 
the main features of the data. The measured kinematic 
dependences will put stringent constraints on any DVCS 
model, and in particular on the generalized parton dis- 
tributions in the nucleon. 
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